Boundary effects in the stepwise structure of the Lyapunov spectra for 
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Boundary effects in the stepwise structure of the Lyapunov spectra and the corresponding wave- 
like structure of the Lyapunov vectors are discussed numerically in quasi-one-dimensional systems 
consisting of many hard-disks. Four kinds of boundary conditions constructed by combinations of 
periodic boundary conditions and hard-wall boundary conditions are considered, and lead to dif- 
ferent stepwise structures of the Lyapunov spectra in each case. We show that a spatial wavelike 
structure with a time-oscillation appears in the spatial part of the Lyapunov vectors divided by mo- 
menta in some steps of the Lyapunov spectra, while a rather stationary wavelike structure appears 
in the purely spatial part of the Lyapunov vectors corresponding to the other steps. Using these 
two kinds of wavelike structure we categorize the sequence and the kinds of steps of the Lyapunov 
spectra in the four different boundary condition cases. 



Q 

U 
d 



(N 
> 

in 

(N 

o 
o 

(N 

o 



> 

x 



I. INTRODUCTION 

Chaos is one of the essential concepts to justify a sta- 
tistical treatment of deterministic dynamical systems. In 
a chaotic system a small initial error diverges exponen- 
tially, as characterized quantitatively by the Lyapunov 
exponents A„, and it means that it is impossible in prin- 
ciple to predict precisely any quantity of deterministic 
systems at any time and a statistical treatment of the 
systems is required. It is well known that even one- 
particle systems can be chaotic and have some of the 
important statistical properties of equilibrium statistical 
mechanics such as the mixing property, etc. For this rea- 
son many works on the subject of chaos have been done 
in one-particle systems, for example, billiard systems and 
Lorentz gas models ^ . However many-particle effects 
should still play an important role in some statistical as- 
pects, such as the central limit theorem and a justifi- 
cation of thermodynamical reservoirs, etc. Therefore it 
should be interesting to know what happens with the 
combination of a chaotic effect and a many-particle ef- 
fect, in other words, to know which chaotic effects do not 
appear in one-particle systems. 

The Lyapunov spectrum is introduced as the sorted set 
{Ai, Aa, • • • } of Lyapunov exponents satisfying the condi- 
tion Ai > A2 > • • • , and is a quantity to characterize 
the many-particle chaotic dynamics. Recently a step- 
wise structure of the Lyapunov spectrum was reported 
numerically, as one of the many-particle chaotic effects, 
in many-disk systems in which each disk interacts with 
the other disks by hard-core interactions || [|, ||. This 
stepwise structure appears in the region of small abso- 
lute values of Lyapunov exponents. This fact suggests 
that steps in the Lyapunov spectra come from slow and 
macroscopic behavior of the system, because small pos- 



itive Lyapunov exponents should correspond to slow re- 
laxation precesses. This point is partly supported by 
the existence of a global structure in the Lyapunov vec- 
tors, the so called Lyapunov modes, which are wavelike 
structures in the tangent space of each eigenvector of a 
degenerate Lyapunov exponent, that is, for each step- 
wise structure Q |[ [| |7| . The wavelike structure of the 
Lyapunov vectors appears as a function of the particle 
position, so this structure is also important as a relation 
to connect the tangent space with the phase space. Al- 
though the Lyapunov vectors have been the subject of 
some works more than a decade ago (for example, see 
Refs. § |, 00 0, HI), it is remarkable that 
an observation of their global structure in fully chaotic 
many-particle systems has only recently appeared. Ex- 
planations for the stepwise structure of the Lyapunov 
spectra have been attempted using periodic orbit models 
p6| and using the master equation approach jlTj . Other 
approaches to the Lyapunov modes have included using 
a random matrix approach for a one-dimensional model 
|fl8f and using a kinetic theoretical approach (n]] , by con- 
sidering these as the " Goldstone modes" . 

If the stepwise structure of Lyapunov spectra is a re- 
flection of a global behavior of the system, then one 
may ask the question: Does such a structure depend 
on boundary conditions which specify the global struc- 
ture of the system? One of the purposes of this paper 
is to answer this question using some simple examples. 
In this paper we investigate numerically the stepwise 
structure of the Lyapunov spectra in many-hard-disk sys- 
tems of two-dimensional rectangular shape with the four 
kinds of the boundary conditions: [Al] purely periodic 
boundary conditions, [A2] periodic boundary conditions 
in the rr-direction and hard- wall boundary conditions in 
the y-direction, [A3] hard-wall boundary conditions in 
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the x-direction and periodic boundary conditions in the 
y-direction, and [A 4] purely hard- wall boundary condi- 
tions. In all cases we took the y-direction as the narrow 
direction of the rectangular shape of the system and in- 
direction as the longer orthogonal direction. The case 
[Al] means that the global shape of the system is like 
the surface of a doughnut shape, the cases [A 2} and [AS] 
means that the system has the shape of the surface of 
pipe with hard-walls in its ends. Case [A 2] is a long pipe 
with small diameter while case [A3] is a short pipe with 
a large diameter. The case [A 4] means that the system 
has just a rectangular shape surrounded by hard-walls. 
Adopting hard-wall boundary conditions in a particular 
direction destroys the total momentum conservation in 
that direction, so by considering these models we can in- 
vestigate the effects of the total momentum conservation 
on the stepwise structure of the Lyapunov spectra and 
the Lyapunov modes. This can be used to check some 
theoretical approaches to these phenomena like those in 
Refs. |l8|; [l^], in which the total momentum con- 

servation plays an essential role in explaining the step- 
wise structure of the Lyapunov spectra or the Lyapunov 
modes. We obtain different stepwise structures of the 
Lyapunov spectra with each different boundary condi- 
tion case. Especially, we observe a stepwise structure 
of the Lyapunov spectrum even in the purely hard-wall 
boundary case [A4 ] , in which the total momentum is not 
conserved in any direction. 

The second purpose of this paper is to categorize the 
stepwise structure of the Lyapunov spectra according to 
the wavelike structure of the Lyapunov vectors. So far 
wavelike structure of the Lyapunov vectors was reported 
in the Lyapunov vector components as a function of po- 
sitions, for example in the quantity Sq^ as a function of 
the position q X j (the transverse Lyapunov mode) [Q, |[ |(J 
and in the quantity &q X j as a function of the position q X j 
(the longitudinal Lyapunov mode) in which Sq^ 

{5q x n j) is the y-component (x-component) of the spatial 
coordinate part of the Lyapunov vector of the j-th parti- 
cle corresponding to the n-th Lyapunov exponent A„, and 
q x j is the ^-component of the spatial component of the j- 
th particle. These wavelike structures appear in the step- 
wise region of the Lyapunov spectrum. However it is not 
clear whether there is a direct connection between the se- 
quence and the kinds of steps of the Lyapunov spectrum 
and the modes of wavelike structure of the Lyapunov vec- 
tors in the numerical studies, that is, how we can catego- 
rize the steps of the Lyapunov spectrum by the Lyapunov 
modes. In this paper we show another type of wave- 
like structure of the Lyapunov vectors, which appears in 
the quantity bqi) Ivyj with the y-component p y j of the 



momentum coordinate of the j-th particle, and use this 
wavelike structure to categorize the stepwise structure 
of the Lyapunov spectra. The Lyapunov vector compo- 
nents Sq X j , Sqljj and Sq^ /p y j have a common feature 
in which each quantity gives a constant corresponding to 
one of the zero Lyapunov exponents of the system with 
purely periodic boundary conditions, because of the con- 
servation of center of mass or the deterministic nature 
of the orbit. More concretely in this paper we consider 
the two kinds of graphs: [Bl] the quantity 8q^ as a 
function of the position q x j, namely the transverse Lya- 
punov mode, and [B2] the quantity 5q^ /Pyj as a func- 
tion of the position q X j and time (or collision number). 
In the two-dimensional rectangular systems consisting of 
many-hard-disks with periodic boundary conditions (the 
boundary case [Al]) it is known that there are two types 
of steps of the Lyapunov spectra: steps consisting of two 
Lyapunov exponents and steps consisting of four Lya- 
punov exponents |^|, The wavelike structure in the 
graph [Bl] corresponding to the two-point steps of the 
Lyapunov spectrum for such a rectangular system is well 
known. In this paper we show the wavelike structure 
corresponding to the four-point steps of the Lyapunov 
spectrum in the graph [B2]. Besides, we observe time- 
dependent oscillations in the graph [B2], whereas the 
graph [Bl] is rather stationary in time. The wavelike 
structure in the graph [B2] also appears in the rectan- 
gular systems with a hard-wall boundary condition (the 
boundary cases [A2], [A3] and [-4^]), specifically even in 
the case of purely hard-wall boundary conditions in which 
the transverse Lyapunov mode [Bl] does not appear. We 
show that the stepwise structure of the Lyapunov spectra 
in the boundary cases [Al], [A 2], [A3] and [A4] can be 
categorized by the wavelike structures of the graphs [Bl ] 
and [B2]. 

One of the problems which make it difficult to inves- 
tigate the structure of the Lyapunov spectra and the 
Lyapunov modes is that the calculation of a full Lya- 
punov spectra for many-particle system is a very time- 
consuming numerical calculation. Therefore it should be 
important to know how we can calculate the stepwise 
structure of the Lyapunov spectra and the Lyapunov 
modes with as little calculation time as possible. It is 
known that a rectangular system has a wider stepwise 
region of the Lyapunov spectrum than a square system 
j5|. Noting this point, in this paper we concentrate on the 
most strongly rectangular system, namely on the quasi- 
one-dimensional system in which the rectangle is too nar- 
row to allow particles to exchange their positions. As will 
be shown in this paper, the stepwise structure of the Lya- 
punov spectrum for the quasi-one-dimensional system is 
the same as the fully rectangular system which allows ex- 
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change of particle positions and each particle can collide 
with any other particle, and the steps of the Lyapunov 
spectra consist of two-point steps and four-point steps 
in the periodic boundary case [Al]. We also determine 
an appropriate particle density to give clearest stepwise 
structure of the Lyapunov spectrum in the quasi-one- 
dimcnsional system. 

The outline of this paper is as follows. In Sect. || we 
discuss how we can get the stepwise structure of the Lya- 
punov spectra for a small number of particle systems with 
a quasi-one-dimensional shape. The density dependence 
of the Lyapunov spectrum for the quasi-one-dimensional 
system is investigated to get a clearly visible stepwise 



structure of the Lyapunov spectra. In Sect. Ill we con- 
sider the purely periodic boundary case (the boundary 
case [j4J]), and investigate wavelike structures in the 
graphs [Bl] and [B2] of Lyapunov vector components. 
In Sect. [fy| we consider the cases including a hard- wall 
boundary condition (the boundary cases [A 2], [A3] and 
[A4]) with calculations of the graphs [Bl] and [B2], and 
compare the stepwise structures of the Lyapunov spectra 
and the wavelike structures in the graphs [Bl] and [B2] 
in the above four boundary cases. Finally we give our 
conclusion and remar ks in Sect. 0. 



II. QUASI-ONE-DIMENSIONAL SYSTEMS 
AND DENSITY DEPENDENCE OF THE 
LYAPUNOV SPECTRUM 



The stepwise structure of the Lyapunov spectra is 
purely a many-particle effect of the chaotic dynamics, 
and so far it has been investigated in systems of 100 
or more particles. However the numerical calculation of 
Lyapunov spectra for such large systems is not so easy 
even at present. Noting this point, in this section we dis- 
cuss briefly how we can investigate the stepwise structure 
of the Lyapunov spectrum for a system whose number of 
particles is as small as possible. We also discuss what is a 
proper particle density to get clear steps of the Lyapunov 
spectrum. These discussions give reasons to choose some 
values of system parameters used in the following sec- 
tions. 
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We consider two-dimensional systems consisting of N 
number of hard-disks in which the radius of the particle 
is R and the width (height) of the system is L x (L y ). One 
way to get the stepwise structure of Lyapunov spectrum 
in a two-dimensional system consisting of a small num- 
ber of particles is to choose a rectangular system rather 
than a square system, because the stepwise region in the 
Lyapunov spectrum is wider in a more rectangular sys- 
tem rather than in a square system even if the number 
of particles is the same p|. Noting this characteristic, 
we concentrate on the most rectangular case, namely the 
quasi-one-dimensional system defined by the conditions 



2RN < L x and 2R < L y < AR. 



(1) 



FIG. 1: Quasi-one-dimensional system: A narrow rectangular 
system satisfying the conditions 2RN < L x and 27? < L y < 
AR. 



A schematic illustration of the quasi-one-dimensional sys- 
tem is shown in Fig. [l| In other words, the quasi-one- 
dimensional system is introduced as a narrow rectangular 
system where each particle can interact with two nearest- 
neighbor particles only and particles do not exchange the 
order of their positions. In the quasi-one-dimcnsional 
system the upper bound p m ax of the particle density p = 
NttR 2 I L x L y is given by p ma x = tt/4 = 0.7853 ••• . In 
such a system we can get a stepwise structure of the Lya- 
punov spectrum even in a 10-particle system (N = 10), 
as shown in Fig. ^|, which is the Lyapunov spectrum nor- 
malized by the maximum Lyapunov exponent Ai « 3.51. 
To get this figure we chose the parameters as R = 1, 
L y = 2Rx(l+ 10~ 6 ), L x = NL V x (1 + 10~ 3 ), and the 
mass M of the particle and the total energy E are given 
by 1 and N, respectively, and we used periodic boundary 
conditions in both the directions. The particle density of 
this system is given by p — 0.7846 • • • . Noting the pairing 
property of the Lyapunov spectrum for Hamiltonian sys- 
tems, namely the property that in Hamiltonian systems 
any positive Lyapunov exponent accompanies a negative 
Lyapunov exponent with the same absolute value fiL , 
we plotted the first half of the Lyapunov spectrum in Fig. 
||. (The same omission of the negative branch of Lya- 
punov spectra from plots will be used throughout this 
paper.) It is clear that the Lyapunov exponents Ai6 and 
A17 form a two-point step in this Lyapunov spectrum. 

In order to calculate Lyapunov spectra and Lyapunov 
vectors we use the algorithm due to Benettin et al., which 
is characterized by intermittent rescaling and renormal- 
ization of Lyapunov vectors p3| , |24| . In the applica- 
tion of this algorithm to systems with hard-core parti- 
cle interactions, we calculate the matrix A(ifc), whose 
column vectors give the Lyapunov vectors 6T n (tk) cor- 
responding to the local-time Lyapunov exponent \ n (tk) 
at time t = tk just after the fc-th collision of parti- 
cles. The dynamics of the matrix A(£&) is given by 
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FIG. 2: Lyapunov spectrum normalized by the maximum 
Lyapunov exponent in a 10-hard-disk system with the peri- 
odic boundary conditions in both directions. The Lyapunov 
exponents Ai6 and A17 shown as the black circles form a 
two-point step accompanying the transverse Lyapunov modes 
shown in Fig. |^. The gray circle is the zero-Lyapunov expo- 
nent A19 whose corresponding Lyapunov vector components 
Sqff show a constant behavior (also see Fig. 0). 



A(tfe-i-i) = NkGkMk^{tk), in which A4 is the matrix 
required to normalize each column vector of the oper- 
ated matrix, Qk is the Gram-Schmidt procedure ensuring 
the orthogonality of the columns of the operated matrix, 
and My°> specifies the tangent space dynamics includ- 
ing a free flight and a particle collision B. The local- 
time Lyapunov exponent \ n (tk) is calculated as the ex- 
ponents of the exponential divergence (contraction) from 
the n-th column vectors of the matrix A(tfc) to the n- 
th column vectors of the matrix GkM-kh-{tk)- The Lya- 
punov exponent A„ is given as a time-averaged local- 
time Lyapunov exponent after a long time calculation: 
A„ = linifc^oo X n (tk). Here we use the standard met- 



ric d6 S * = ELldSqS 2 + d6q { y f 2 + d6p^ 2 + dS Py f 2 } 
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and 



for the tangent space with the x-component 5p XJ 
y-component Sp y n ^ of the Lyapunov vector of the j-th 
particle corresponding to the local-time Lyapunov expo- 
nent A„. Other articles such as Refs. J3|, [2^, Q should be 
referred to for more details of the Benettin algorithm and 
the tangent space dynamics of many-hard-disk systems. 

It is very important to note that there are two ways 
that convince us of the structure of the Lyapunov spec- 
trum; one is simply to find a step structure directly in the 
Lyapunov spectrum, and another is to find a structure 
in the Lyapunov vectors corresponding to specific Lya- 
punov exponents. Fig. ^ is the plot of the time-average 



(Sq^j) of the y-component of the j-th particle of the 
spatial part of the Lyapunov vector ST n corresponding 
to the Lyapunov exponents A„, for n = 19, 17 and 16 as 
functions of the time- average (q X j)/L x of the normalized 
x-component of spatial coordinate of the j-th particle 
(j = 1,2, ••• ,N) (the graph [Bl]). The corresponding 
Lyapunov exponents A17 and Ai6 are shown as the black- 
filled circles in Fig. |^ and construct a two-point step of 
the Lyapunov spectrum. To take the time average of the 
quantities 5q y n ^ and q X j we picked up these values just 
after particle collisions over 1007V collisions and take the 
arithmetic averages of them. The step consisting of two 
points in the Lyapunov spectrum of Fig. |^ accompanies 
wavelike structures in their Lyapunov vectors, which are 
called the transverse Lyapunov modes. It should be em- 
phasized that the Lyapunov modes in Fig. [| are rather 
stationary over 1007V particle collisions and to take the 
time average can be a useful way to get clear their wave- 
like structures. In this figure we fitted the numerical 
data for the Lyapunov modes corresponding to the Lya- 
punov exponents A n , n — 17 and 16, to the function 
y = a n cos(27re + (3 n ) ((a 17 ,/3i 7 ) = (0.19913,0.9703) for 
the triangle dots, and (a 16l Pi 6 ) = (0.24025,-0.59677) 
for the square dots). It should be noted that the differ- 
ence /?i7-/3i 6 = 0.9703-(-0.59677) = 1.56707 of the two 
values of the phases (3 n is approximately 7r/2 = 1.570 • ■ ■ , 
meaning that these two waves are orthogonal with each 
other. The graph of the Lyapunov vector component 
5q y n ^ as a function of the position q x j (j = 1, 2, • • • , TV) 
become constant in one of the zero Lyapunov exponent 
A19 shown as the gray-filled circle in Fig. g. In Fig. || we 
also fitted the numerical data corresponding to the zero- 
Lyapunov exponent A19 by the constant function y = aig 
with the fitting parameter value a±g = 0.015114. 



FIG. 3: Time-averaged j/-components (Sq^ ) of the spatial 
part of the Lyapunov vector of the j-th particle as functions 
of the time-average {q x j)/L x of the normalized z-component 
of spatial coordinate of the j-th particle corresponding to the 
Lyapunov exponents X n , n = 19, 17 and 16, in the 10-hard- 
disk system with the periodic boundary conditions in both 
directions. The circle, triangle and square dots correspond to 
the Lyapunov exponents A19, A17 and Ai6, respectively .which 
are shown as the black- and gray-filled circles in Fig. PJ. The 
dotted and broken lines are the fitting lines for the sinusoidal 
functions, and the solid line is the fitting for the constant 
function. 



Although we can recognize the two-point step of the 
Lyapunov spectrum in Fig. |], it is important to note 
that there is another type of steps in the two-dimensional 
hard-disk system with a rectangular shape and periodic 
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boundary conditions. It is well known that in such a sys- 
tem the Lyapunov spectrum can have a stepwise struc- 
ture consisting of two-point steps and four-point steps 
H 19 1 . If we want to investigate the four-point steps 
of the Lyapunov spectrum we may have to consider a 
system consisting of more than 10 particles. 

Now we mention an angle in components of the Lya- 
punov vectors. Fig. |] is the graph for the time-average 
(6 n ) of the angle 9 n between the spatial part <5q(") 
and the momentum part Sp^ of the Lyapunov vector 
5T n = (5q^ n \ Sp^) T where T is the transverse opera- 
tion, and n is defined by 



cos 



V|(5q(")||<Sp(") 



(2) 



except for the ones corresponding to the zero Lyapunov 
exponents. Here in order to take the time-average of 
the angle 6 n we picked up its values just after particle 
collisions over 10007V collisions and took the arithmetic 
average of them. This graph shows that the spatial part 
<5q(") and the momentum part 5p^ of the Lyapunov 
vector are pointed toward almost the same direction, es- 
pecially for the Lyapunov vectors corresponding to small 
absolute values of Lyapunov exponents This fact 

suggests that if we get a structure in the vector Sq^ 
then we may get the similar structure in the vector Sp^ , 
and vice versa. It should also be noted that this gives a 
justification for some approaches to Lyapunov exponents 
in which the Lyapunov exponents are calculated through 
the spatial coordinate part only (or the momentum part 
only) of the Lyapunov vector |H| |25|] . We may also note 
that the graph for 6 n , n = 1, 2, • • ■ , 2N — 3 and ^4jv-n+i, 
n = 1, 2, • ■ • , 2N — 3 are symmetric with respect to the 
line y = 1/2, which comes from the symplectic property 
of the Hamiltonian dynamics. 



FIG. 4: Time-averaged angle (#„) between the spatial part 
<5q'"' and the momentum part Jp 1 -™' of the Lyapunov vector 
ST n = (<5q' n \ <5p' n ') T corresponding to the Lyapunov expo- 
nent A„ in the 10-hard-disk system with the periodic bound- 
ary conditions in both directions. 



Another important point to get a clear stepwise struc- 
ture for the Lyapunov spectrum in a small system is to 
choose a proper particle density of the system. Even 
if we restrict our consideration to quasi-one-dimensional 
systems, the shape of the Lyapunov spectrum depends on 
the particle density, so we should choose a density that 
gives a clearly visible stepwise structure in the Lyapunov 
spectrum. Fig. |^(a) is the Lyapunov spectra normal- 



ized by the maximum Lyapunov exponents for quasi-one- 
dimensional systems consisting of 50 hard-disks (N = 50) 
with periodic boundary conditions in both directions. We 
also give an enlarged figure |^(b) for the small Lyapunov 
exponent region. Here the system parameters are given 
by R = 1, L y = 2i?(l + 10~ 6 ), L x = NL y (l + d), and we 
used the case of M = 1 and E = N. The five Lyapunov 
spectra correspond to the case of d = 10~ 4 (the circle 
dots, the density p — 0.7853 • • •), the case of d = 10~ 4 
(the triangle dots, the density p — 0.7139 • • • ), the case 
of d = 1 (the square dots, the density p = 0.3926 • • ■ ), 
the case of d = 10 (the diamond dots, the density 
p = 0.07139 ■ • • ), the case of d = 10 2 (the upside-down 
triangle dots, the density p = 0.007776 • ■ •). The max- 
imum Lyapunov exponents Ai are given approximately 
by 3.62, 2.44, 0.934, 0.279 and 0.0579 in the cases of 
d = 10~ 4 , 10 _1 , 1, 10 and 10 2 , respectively. As shown 
in Figs. |^(a) and (b), for smaller values of the quan- 
tity d (namely in the higher particle density) the gaps 
between the nearest steps of the Lyapunov spectrum be- 
come larger, although the stepwise region of the Lya- 
punov spectrum does not seem to depend on the quantity 
d. This means that we can get a clear stepwise structure 
of the Lyapunov spectrum in the small d case (namely at 
high density). 



FIG. 5: Density dependence of Lyapunov spectrum nor- 
malized by the maximum Lyapunov exponent for quasi-one- 
dimensional systems consisting of 50 particles with the peri- 
odic boundary conditions in both directions. The five Lya- 
punov spectra correspond to the case of L y = 2R(l + 10 -6 ), 
L x = NL y (l + d) with d = 10~ 4 (the circle dots), the case 
of d = 1CP 1 (the triangle dots), the case of d — 1 (the square 
dots), the case of d = 10 (the diamond dots), the case of 
d — 10 2 (the upside-down triangle dots), (a) Full scale, (b) A 
small positive Lyapunov exponent region including the step- 
wise structure of the Lyapunov spectra. In the figure (b) 
we filled dots by black (gray) for the Lyapunov exponents 
whose corresponding time-averaged Lyapunov vector compo- 
nents {Sq^ ) of the spatial part of the Lyapunov vector of the 
j-th particle as functions of the time-averaged a;-component 
{q X j) of spatial coordinate of the j-th particle show wavelike 
structures (constant behaviors). 



In Fig. ||(b) the Lyapunov exponents accompany- 
ing wavelike structures (constant behaviors) in the time- 
averaged Lyapunov vector components (^Q^)i as func- 
tions of the position (q X j), are shown as the dots blacked 
(grayed) out. In the small d case, looking from the zero 
Lyapunov exponents, the two-point step appears first 
(see the cases of d = 10~ 4 , 10~ 4 and 1 in Fig. ||(b)), 
whereas the four-point step appears first in the large d 
case (see the cases of d = 10 1 and d = 10 2 in Fig. ||(b)). 
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Besides, at least in the small d case, the two-point steps 
and the four-point steps do not appear repeatedly (see 
the cases of d = 1CT 4 and 1CT 1 in Fig. |(b)). These 
facts mean that the sequence of steps in the Lyapunov 
spectrum depends on the quantity d, namely the particle 
density. 



Wavelike structures in the time-averaged Lyapunov 
vector components (5q y j) as functions of the position 
(q x .j), namely the transverse Lyapunov modes, appear 
mainly in two-point steps of the Lyapunov spectra. 
Therefore we can use these wavelike structures to distin- 
guish two-point steps from four-point steps of the Lya- 
punov spectra. However such a distinguishability some- 
times seems to fail in steps of the Lyapunov spectra 
near a region where the Lyapunov spectra is changing 
smoothly. Actually the transverse Lyapunov modes may 
appear even in some four-point steps, if they are near 
such a smoothly changing Lyapunov spectrum region. 
On the other hand, in such a region of the Lyapunov 
spectra, fluctuation of Lyapunov vectors is rather large, 
and the wavelike structure of the Lyapunov vectors be- 
come vague. In Fig. ||(b) we did not indicate by black- 
filled dots the Lyapunov exponents whose corresponding 
Lyapunov vector components {^QZj ) &s functions of the 
position (q x j) show wavelike-structures with such a level 
of vagueness. 



Based on the discussions in this section, in the follow- 
ing two sections we consider only the case whose system 
parameters are given by N = 75, R = 1, M = 1 and 
E = N . The height and the width of the system are 
given by L y = 2R(1 + 1(T 6 ) and L x = 1.5NL y (the den- 
sity p — 0.5235 • • • ) in the purely periodic boundary case 
considered in the next section. In this will be 

shown in the next section, we can recognize at least 2 
clear four-point steps, and the two-point steps and the 
four-point steps of the Lyapunov spectrum appear re- 
peatedly in the first few steps in the Lyapunov spectrum 
for the system with the purely periodic boundary con- 
ditions. We always take the time averaged quantities 
) an< ^ (fei) °f t ne quantities Sq^ and q X j, respec- 
tively as the arithmetic average of their values taken in 
the times immediately after particle collisions, over 1001V 
collisions. We calculated more than 2 x 10 5 collisions (10 6 
collisions in some of the models) in order to get the Lya- 
punov spectra and the Lyapunov vectors in the models 
considered in this paper. 



III. QUASI-ONE-DIMENSIONAL SYSTEMS 
WITH PERIODIC BOUNDARY CONDITIONS 

In this section we consider the Lyapunov spectrum for 
the quasi-one-dimensional system with periodic bound- 
ary conditions in both the directions (the boundary case 
[j4J]). A schematic illustration of this system for latter 
comparisons is given in Fig. |^ in which the broken line of 
the boundary means to take the periodic boundary con- 
ditions. This system satisfies total momentum conserva- 
tion in both the directions, and is regarded as a reference 
model for the models considered in the following section. 

iO"0"0"f"0"0: 

1) 

FIG. 6: A schematic illustration of a quasi-one-dimensional 
system with periodic boundary conditions in both directions. 
The broken lines indicate periodic boundary conditions at 
that boundary. 

Fig. [?] is the small positive Lyapunov exponent region 
of the Lyapunov spectrum normalized by the maximum 
Lyapunov exponent Ai w 1.33, including its stepwise re- 
gion, for a quasi-one-dimensional system with periodic 
boundary conditions in both the directions. The global 
shape of the Lyapunov spectrum is given in the inset in 
this figure. We used the values of the system parameters 
chosen at the end of the preceding section. At least 5 
steps consisting 3 two-point steps and 2 four-point steps 
are clearly visible in this Lyapunov spectrum. 

FIG. 7: The stepwise structure of the Lyapunov spectrum 
normalized by the maximum Lyapunov exponent for a quasi- 
one- dimensional system with periodic boundary conditions in 
both directions. Inset: Full scale for the normalized Lya- 
punov spectrum. The circular dots are filled by black (gray) 
in the Lyapunov exponents corresponding to wavelike struc- 
tures (a constant behavior) of the time-averaged Lyapunov 
vector components (5<j£™ ) shown in Fig. |[ 

The two-point steps of the Lyapunov spectrum ac- 
company wavelike structures in their corresponding Lya- 
punov vectors. Fig. ^ is the graphs of the time-averaged 
Lyapunov vector components {Sq y j} corresponding to 
the Lyapunov exponent A n , n = 149,147,146,141 and 
140 as functions of the time-averaged position compo- 
nent (q X j) normalized by the length L x (the graph [-81]). 
The Lyapunov exponents used for this figure are shown 
as the black or gray circles in Fig. (?]. In this figure 
we also gave fittings of the numerical data by sinusoidal 
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equations and a constant function. The fitting equations 
are y = a n cos(27ra; + /? n ) for n = 147, 146 ( (01147, /?i47) = 
(-0.16132,-4.5575), (a U 6, /?i4e) = (-0.16125,0.14928) 
), and y = a n cos(47rx + (3 n ) for n = 141,140 ( 
(oi4i,/5i4i) = (-0.15775,-0.30581), (o 140 ,/5i4o) - 
(0.15782,1.2666) ). Concerning these fitting parame- 
ters, it is important to note the relations |oi47| ~ |«i46|; 
oi4i| ~ |«i4o|, /?i46-/3i47 ~ 3tt/2 and fluo - Pm « 7r/2, 
meaning that the two wavelike structures in the same 
two-point step are orthogonal to each other. In Fig. || 
we also draw a graph of the averaged Lyapunov vector 
component as a function of the normalized posi- 

tion {q X j)/L x corresponding to the second zero Lyapunov 
exponent A149, which is roughly speaking constant, and is 
fitted by a constant function y = 0149 with the fitting pa- 
rameter value ai49 = 0.0 1 8901. These results suggest a 
conjecture for the approximate form of the Lyapunov vec- 
tor component Sq^ k ^ corresponding to the Lyapunov 
exponents \ n ( k ) and \i(k)-x m the same fc-th two-point 
step counting from the zero Lyapunov exponents as 

f /27Tfc \ . /27Tfc \\ 

< a k cos I -j—q X j + Pk J , a k sm I ——q X j + Pk J > , 

(3) 

j = 1,2, ••• ,N with constants a k and f3 k - (Note that 
in this paper we always count the sequence of steps of 
Lyapunov spectra looking from the zero Lyapunov expo- 
nents, so, for example, the two-point step consisting of 
the Lyapunov exponents A147 and A146 (A141 and A 140) in 
the Lyapunov spectrum shown in Fig. || is the first (sec- 
ond) two-point step.) It may be noted that the constant 
Ofc in Eq. (|3|) should be determined by the normalization 
condition of the Lyapunov vector, which is used in the 
Benettin algorithm to calculate the Lyapunov spectrum 
and Lyapunov vectors in this paper. 

FIG. 8: Time-averaged Lyapunov vector components {5<& ) 
as functions of the time-average (q X j)/L x of the normalized 
x-component of the spatial coordinate of the j-th particle cor- 
responding to the n-th Lyapunov exponent X n for the quasi- 
one-dimensional system with the periodic boundary condi- 
tions in both directions (n = 149, 147, 146, 141 and 140). The 
corresponding Lyapunov exponents are shown as the black- 
and gray-filled circles in Fig. ^. The numerical data are fitted 
by a constant function and sinusoidal functions. 

It should be emphasized that the wavelike structures 
in Figs. H and || correspond to the two-point steps of 



the Lyapunov spectrum. We cannot recognize such a 
clear wavelike structure in the graph of the Lyapunov 
vector component (dq^) as a function of the position 
(q x j) corresponding to the four-point steps of the Lya- 
punov spectrum. This suggests that the physical mean- 
ing of the four-point steps is different from the two-point 
steps, and we should consider a different quantity to in- 
vestigate as a structure of Lyapunov vectors correspond- 
ing to the four-point steps of the Lyapunov spectrum. 
Now, as one of the important results of this paper, we 
show that wavelike structures in the quantities Sq y r ^ /p y j 
as functions of q X j and the collision number appear in 
the four-point steps of the Lyapunov spectrum. As a 
motivation to introduce the quantity $q y j /p y j to inves- 
tigate its structure, we note that the small perturbation 
of the spatial coordinates of the particles in the direc- 
tion of the orbit, namely the small perturbation 5q oc p, 
leads to a zero- Lyapunov exponent, so that the quantities 
Sqyj /Pyj, j = 1, 2, • ■ • , N should give a constant value 
corresponding to this zero Lyapunov exponent. This is 
the common feature as the quantity 8q~"j , which shows 
a constant behavior corresponding to one of the zero- 
Lyapunov exponents owing to the conservation of center 
of mass, and whose wavelike structure we have already 
discussed. 

Before showing graphs of the quantities dq^' /p y j, we 
discuss some difficulties in the investigation of wavelike 
structure in these quantities. It is much harder to get the 
wavelike structure of these quantities corresponding to 
the four-point steps of the Lyapunov spectrum compared 
to the quantities 5q^ corresponding to the two-point 
steps of the Lyapunov spectrum, for at least two reasons. 
As the first reason, fluctuation in the wavelike structure 
of the quantities 5q y j /p y j is much bigger than in the 
wavelike structure of the quantities Sq^ , partly because 
such a fluctuation is magnified in the case of a small ab- 
solute value of the quantity p y j appearing in the denom- 
inator of the quantity Sq y ^ I Pyj- Secondly, the wavelike 

structure of the quantity Sqlff /Pyj, at least the magni- 
tude of its wave, oscillates periodically in time, whereas 
the wavelike structure of the quantities Sq y ^ in Fig. || 
are stationary at least over more than 100 A*" collisions. 
This fact gives an upper bound on the time period (or 
the collision number interval) over which we can take the 
time-average of the quantities ^qff I Pyj in order to sup- 
press the fluctuations and get their clear wavelike struc- 
tures. In this paper we express the local time-averages 
of the quantities Sqy/Pyj and q X j as (Sq y j /Pyj)t and 
(q x j)t, respectively, with the suffix t to remind us that 
they can change in time. 

In this and the next sections we will give the many 
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FIG. 9: Local time-averaged quantities (<5<7y™ /Pvj)t as func- 
tions of the normalized position (q x j)t/L x and the collision 
number n t corresponding to (a) the Lyapunov exponent A142 
in the first four-point step and (b) the Lyapunov exponent 
A139 in the second four-point step in the same collision num- 
ber interval [543000,548100]. The system is the quasi-one- 
dimensional system with the periodic boundary conditions in 
both directions, and the corresponding Lyapunov exponents 
are indicated by arrows in Fig. M. In the contour plots on the 
bottoms of these three-dimensional plots, red dotted lines, 
black solid lines and blue broken lines correspond to the val- 
ues {Sq^jf /pyj)t = 0.08,0 and —0.08, respectively. 



graphs of <5<zi" I Pyj as functions of q X j and the collision 
number by taking their local time-averages, so here we 
summarize how we calculate the data for those graphs 
from a technical point of view. First we take the arith- 
metic time-average (Sq y j /pyj)t and (q X j)t of the quan- 
tities Sq^/p y j and q x j, respectively, using their values 



just after particle collisions over 47V collisions (8N col- 
lisions), but if the absolute value \p y j\ of the momen- 
tum is less than 10 percent (5 percent) of the averaged 
momentum amplitudes y/2ME/N then we exclude the 
quantity Sq^ /p y j at that time from samples to take 
this local time-average, in the models of this section and 
the subsection pV A (in the models of the subsections 
1VB and IV Cj ). (Therefore the sample number for tak- 
ing the arithmetic averages can be less than AN (8N) 
in the models of this section and the subsection |IV A 
(in the models of the subsections IV B and JVC).) On 
this local time-averages we can still get more than ten 
locally time-averaged datas for different times in one pe- 
riod of the time-oscillations for the graphs in the slow- 
est time-oscillating movement of the quantity 5q y n ^ /p y j 
corresponding to a step of the Lyapunov spectra, for ex- 
ample, corresponding to the first four-point step of the 
model in this section. In this paper we consider the graph 
of the quantity {Sq y ^ /Pyj)t as a function of {q X j)t and n t 
which is the first collision number of the collision number 
interval taking the local time averages (•••)*■ 

Figure ^(a) and (b) is the graphs of the quantity 
(Sq y ^ /p y j)t as functions of the normalized position 
{q X j)t/L x and the collision number n t (the graph [B2]), 
corresponding to the Lyapunov exponents X n , n = 145 
and 139, respectively, indicated by arrows in Fig. |t]. 
These two graphs correspond to Lyapunov exponents 
in different four-point steps and are given in the same 
collision number interval [543000,548100]. In the graph 
corresponding to the Lyapunov exponent A145 (A139) we 
can recognize a spatial wavelike structure with the spa- 
tial wavelength 1 (1/2) oscillating in time. The time- 



oscillating period corresponding to the Lyapunov ex- 
ponent A139 is about half of the period of the oscil- 
lation corresponding to the Lyapunov exponent A145. 
These graphs, especially Fig. |](b), are the most diffi- 
cult graphs to recognize the structures from their three- 
dimensional plots in this paper, and in order to recog- 
nize the structures the contour plots given in the bot- 
toms of these three-dimensional plots may be helpful. 
In these contour plots we color-coded mountain regions 
((Sq y f /p y j)t > 0) and valley regions ((Sq y f /p vj ) t > 0) 
of the time-oscillating wavelike structure of the quan- 
tities (Sq y ^ /p y j)t by the red dotted contour lines of 
(filyf /Pyj)t = 0.08 and the blue broken contour lines 
of {Sq y " /p y j)t — —0.08, respectively, separating them 
by the black solid contour lines of {5q^f /p y j)t — 0. 



Now we consider a relation among the quantities 
(5q y ^ /p V j)t corresponding to Lyapunov exponents in the 



same four-point step. Fig. 10 is the contour plots of the 
quantity (5q y j /p y j)t as functions of the normalized po- 
sition {q X j)t/L x and the collision number n t correspond- 
ing to the first four-point step consisting of the Lyapunov 
exponents A„, n — 145, 144, 143 and 142 in the same col- 
lision number interval [538500,548100]. Here red dotted 
lines, black solid lines and blue broken lines correspond 
to (Sq y j /p y j)t = 0.08,0 and —0.08, respectively. It is 
clear that in these four graphs their spatial wavelengths 
(given by the system length L x ) and time-oscillating pe- 
riods (given by a constant 7o) are almost the same as 
each other. On the other hand we can recognize that 
the position of the nodes of the spatial waves of the 
graphs 10(a) and |l(](b), as well as the nodes of the spa- 
tial waves of the graphs 0(c) and |o|(d), coincide with 
each other approximately, and the phase difference be- 
tween the spatial waves of the graphs |l^(a) and |To|(c) 
is 7r/2 approximately. Besides, the node of the time- 
oscillation of the wave-amplitude of the graphs [Tc|(a) 
and |l0|(c), as well as the nodes of the time-oscillation 
of the wave-amplitude of the graphs |l^(b) and |Toj(d) , 
coincide with each other approximately, and the phase 
difference between time-oscillations of the graphs |lC|(a) 
lC(b) is about 7r/2. These points are summarized in Fig. 
1 1 , which is the schematic illustration of the phase re- 
lations among the time-oscillating wavelike structures of 
the quantity (5q y f /p yj ) t , Here the phase [PI], [P2], [PS] 
and [P4] correspond to Figs. |l0](a), (b), (c) and (d), 
respectively. These observations suggest that the Lya- 

4. j- X X (™( fe )-!) X (™( fe )- 2 ) 

punov vector components oq y j , oq y j , aq y j 

and Sq y j ^ corresponding to the Lyapunov exponents 
constructing the fc-th four-point step are approximately 
expressed as 
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FIG. 10: Contour plots of the local-time averaged quantities {Sq y " /pyj)t as functions of the normalized position (q X j)t/L x and 
the collision number n t corresponding to the first four-point step consisting of the Lyapunov exponents A n , n — 145, 144, 143 
and 142 in the same collision number interval [538500, 548100]. Here red dotted lines, black solid lines and blue broken lines 
correspond to the values {Sq y j /p V j)t = 0.08,0 and —0.08, respectively. The system is the quasi-one-dimensional system with 
the periodic boundary conditions in both directions, and the corresponding Lyapunov exponents used in this figure are indicated 
by the brace under circles in Fig. R. 



fxJHk)) . (n(fe)-l) c (n(fc)-2) ,(h{k)-3)\ 

\ uq yj ' uq yj ' uq yj ' uq yj j 



f2nk - \ f2nk . . 
otkPyjCOS \-j—q X j+Pk I cos \ -^nt+jk ) , 

/2irk ~ \ ( 2-nk 
otkPyjCOS [-j—'lrj I >';. j sin ( ^— j 

/ 2irk - \ ( 2ixk 

sin 



& k pyjsm y-j—lxj+Pkj cos ) , 

sin ( -j— q X j+(3k) sin I -^n t +^ k 



(4) 



j = 1,2,- •• ,7V, with constants ak, Pk and 7fc. Here 
7q is the period of the time-oscillation of the quantity 
^Qyj ^ /Pyj m ^ ne nrs ^ four-point step. (In this paper we 
use the quantity 7o as the period of the particle-particle 
collision number, but we can always convert it into the 
real time interval approximately by multiplying it with 
the mean free time, which is, for example, about 0.0248 
in the model of this section.) It may be emphasized that 
the level of approximation in Eq. (4) for the four-point 
steps may be worse than in Eq. (3) for the two-point 
steps. 



IV. QUASI-ONE-DIMENSIONAL SYSTEMS 
WITH A HARD- WALL BOUNDARY CONDITION 



+ <&<!,/■> /P,j>, > 
_ < 59v/"»/Pv/>, = 
- <M,!->/P,i>i < 



Phase [ PI ] 



Phase [ P2 



Phase [ P3 



Phase [ P4 



n, n, + To n, + 2To n, n, + To n, + 2To 

FIG. 11: Schematic illustration of the phase relations among 
the four time-oscillating wavelike structures of the quantities 
(Sq^ /Pyj)t as functions of the position (q x j)t and the colli- 
sion number nt, corresponding to the first four-point step in 
the same collision number interval. The spatial wavelength 
and the period of the time-oscillations are given by L x and 
To, respectively. Thick gray lines means node lines, the region 
indicated by a red plus sign (+) is the region where the quan- 
tity {&Qyf /Pyj)t is positive, and the region indicated by a blue 
minus sign (— ) is the region where the quantity (Sq^ /p y j)t 
is negative. The phases [PS] and [P4] differ from the phases 
[PI] and [P2] by a phase shift of 7r/2 in the spatial direction. 
The phases [PI] and [P3] differ from the phases [P2] and [P4] 
by a phase shift of n/2 in time. 



In this section we consider quasi-one-dimensional sys- 
tems with boundary conditions including a hard-wall 
boundary condition. Noting that there are two directions 
in which we have to introduce the boundary conditions, 
we consider the three cases of systems; [A 2] the case of 
periodic boundary conditions in the cc-direction and hard- 
wall boundary conditions in the y-direction, [A3] the case 
of hard-wall boundary conditions in the x-direction and 
periodic boundary conditions in the y-direction, and [A4] 
the case of hard-wall boundary conditions in both the 
directions. To adopt the hard-wall boundary condition 
means that there is no momentum conservation in the 
direction connected by hard-walls, so it allows us to dis- 
cuss the role of momentum conservation in the stepwise 



structure of the Lyapunov spectrum by comparing mod- 
els having a hard-wall boundary conditions with models 
having periodic boundary conditions. We will get differ- 
ent stepwise structures for the Lyapunov spectra in the 
above three cases and the case of the preceding section, 
and the investigation of Lyapunov vectors corresponding 
to steps of the Lyapunov spectra leads to relating and to 
categorizing them with each other. 

In the systems with a hard-wall boundary condition we 
should carefully choose the width L x and the height L y of 
the systems for meaningful comparisons between the re- 
sults in the systems with different boundary conditions. 
It should be noted that in systems with pure periodic 
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boundary conditions, the centers of particles can reach 
to the periodic boundaries, while in hard-wall boundary 
conditions the centers of particles can only reach within a 
distance R (the particle radius) of the hard-wall bound- 
aries. In this sense the effective region for particles to 
move in the system with hard-wall boundary conditions 
is smaller than in the corresponding system with peri- 
odic boundary conditions, if we choose the same lengths 
L x and L y . In this section the lengths L x and L y of the 
systems with a hard-wall boundary condition are cho- 
sen so that the effective region for particle to move is 
the same as in the pure periodic boundary case consid- 
ered in the previous section, and are given by (L y , L x ) = 
(2R(l + 10~ 6 ) + 2R,1.5N(L y - 2R)) in the case [A2], 
(L y , L x ) = (2R(l + 10~ 6 ), l.hNLy + 2R) in the case [A3] 
and (L y , L x ) = (2R(l + 10- 6 ) + 2R,1.5N(L y -2R) + 2R) 
in the case [A4 ] . These choices of the lengths L x and L y 
also lead to the almost same mean free time in the four 
different boundary condition cases. In the cases of [A2] 
and [A4] with this choice of the system width L y , in prin- 
ciple there is a space to exchange the particle positions 
in the x-direction (namely, strictly speaking these cases 
do not satisfy the second condition of Eq. (0) ), but the 
space is extremely narrow (that is 2R x 10" 6 ) so that it 
is almost impossible for particle positions to actually be 
exchanged. 

A. The case of periodic boundary conditions in the 
T-direction and hard-wall boundary conditions in 
the y-direction 

The first case is the quasi-one-dimensional system with 
periodic boundary conditions in the ^-direction and with 
hard-wall boundary conditions in the y-direction (the 
boundary case [^4£])- A schematic illustration of this 
system is given in Fig. ^ in which periodic boundary 
conditions and hard-wall boundary conditions are repre- 
sented as bold solid lines and broken lines, respectively. 

• • o 'I; Uti 

FIG. 12: A schematic illustration of a quasi-one-dimensional 
system with periodic boundary conditions in the x-direction 
and hard-wall boundary conditions in the {/-direction. The 
broken line on the boundary means to take periodic boundary 
conditions, and the solid line on the boundary means to take 
hard-wall boundary conditions. 

Figure |l3| is the Lyapunov spectrum normalized by the 
maximum Lyapunov exponent Ai ss 1.30 for this sys- 



FIG. 13: Stepwise structure of the Lyapunov spectrum nor- 
malized by the maximum Lyapunov exponent for the quasi- 
one- dimensional system with periodic boundary conditions in 
the :r-direction and hard-wall boundary conditions in the y- 
direction. Inset: Full scale of the normalized Lyapunov spec- 
trum. 



tern. In this figure we showed a small positive region 
of the Lyapunov spectrum including its stepwise struc- 
tures, while the full scale of the positive branch of the 
Lyapunov spectrum is shown in the inset. In this sys- 
tem the y-component of the total momentum is not con- 
served because of the hard-wall boundary conditions in 
the y-direction, so there are only 4 zero-Lyapunov ex- 
ponents in this system. This figure shows clearly that 
the steps of the Lyapunov spectrum consist of four-point 
steps only, and there is no two-point step in the Lya- 
punov spectrum which appears in the model discussed 
in the previous section. Besides, we cannot recognize a 
wavelike structure in the graph of Lyapunov vector com- 
ponents {Sq y j) as a function of the position {q X j) (the 
graph [Bl]) in this model. A comparison of this fact with 
results in the previous model suggests that the two-point 
step of the Lyapunov spectrum in the previous section 
should be strongly connected to the conservation of the 
y-component of the total momentum. 



FIG. 14: Local time-averaged quantities (5<j£" /Pyj)t as func- 
tions of the normalized position (q x j)t/L x and the collision 
number n t corresponding to the Lyapunov exponents Ai48 and 
A144 in the same collision number interval [390600,396000]. 
The system is the quasi-one-dimensional system with peri- 
odic boundary conditions in the a;-direction and hard-wall 
boundary conditions in the y-direction, and the correspond- 
ing Lyapunov exponents are indicated by arrows in Fig. [b| 
Contour plots on the bottoms of these three-dimensional plots 
are given by red dotted lines, black solid lines and blue bro- 
ken lines corresponding to the values {Sq^ /p y j)t = 0.08,0 
and —0.08, respectively. 

We consider a relation between the four-point steps 
in the model of this section and in the model of the 
previous section by investigating the graph of the quan- 
tity ($C[yj* /Pyj)t as a function of the normalized position 
(q X j)t/L x and the collision number n t (the graph [-B.S]). 
Fig. |l4] is such graphs, and corresponds to the Lya- 
punov exponent A148 in the first four-point step ( Fig. 
|lj(a) ) and the Lyapunov exponent A144 in the second 
four-point step ( Fig. [l4](b) ), which arc indicated by 
arrows in Fig. [l^, in the same collision number interval 
[390600, 396000]. The wavelike structures of these graphs 
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FIG. 15: Contour plots of the time-averaged quantities (Sq^ /Pyj)t as functions of the normalized position (q X j)t/L x and 
the collision number n t corresponding to the first four-point step consisting of the Lyapunov exponents A n , n — 148, 147, 146 
and 145 in the same collision number interval [385800, 396000]. The system is the quasi-one-dimensional system with periodic 
boundary conditions in the a;-direction and hard-wall boundary conditions in the y-direction, and the corresponding Lyapunov 
exponents are indicated by the brace under circles in Fig. Il3l Here red dotted lines, black solid lines and blue broken lines 
correspond to the values {8q£j' /p v j)t = 0.08,0 and —0.08, respectively. 



have a wavelength l/i for the z-th four-point steps. The 
time-oscillating period corresponding to the Lyapunov 
exponent Ai48 is almost the same as the period Tq of the 
first four-point steps of the previous model, and is ap- 
proximately twice as long as the time-oscillating period 
in the Lyapunov exponent A144 of this model. These fea- 
tures are common with the four-point steps in the mod- 
els of the previous section, suggesting that the four-point 
steps of the Lyapunov spectrum in this model correspond 
to the four-point steps in the model of the previous sec- 
tion. We can also show that the quantities (Sq y j /Pyj)t 
corresponding to the zero Lyapunov exponents A 150 and 
A149 are constants as functions of (q X j)t and nt approxi- 
mately. 



B. The case of hard- wall boundary conditions in 
the x-direction and periodic boundary conditions in 
the y-direction 

As the next system we consider a quasi-one- 
dimensional system with hard-wall boundary conditions 
in the ^-direction and periodic boundary conditions in 
the y-direction (the boundary case [.A 5]). A schematic 
illustration of such a system is given in Fig. ^ with solid 
lines for the hard-wall boundary conditions and broken 
lines for the periodic boundary conditions. 



Now we proceed to investigate the graph of the quan- 
tities (Sqffi /Pyj)t corresponding to the Lyapunov ex- 
ponents in the same four-point steps of the Lyapunov 
spectrum. Fig. [l5| is the contour plots of these quan- 
tities as functions of the normalized position (q X j)t/L x 
and the collision number n t , corresponding to the first 
four-point step consisting of the Lyapunov exponents A„, 
n = 148, 147, 146 and 145, in the same collision number 
interval [385800,396000]. The corresponding Lyapunov 
exponents are indicated by the brace under circles in Fig. 
|l3| . In this case the contour lines of (Sqff /p y j)t — 
seem to be slanting, but if we pay attention of the moun- 
tain regions and the valley regions in these graphs then 
we can realize that the structure of these four graphs 
are similar to the contour plots of the four graphs in 
Fig. |l^ for the previous model. Therefore the phase 
relations among the time-oscillating wavelike structures 
of the quantity (Sq!$ /Pyj)t can be summarized in the 
schematic illustration given in Fig. ^ like in the pre- 
vious model. This suggests an approximate expression 
for the Lyapunov vector components given by Eq. (Q), 
for Lyapunov vector components Sq^ corresponding to 
the Lyapunov exponents of the four-point steps in this 
model. 



FIG. 16: A schematic illustration of a quasi-one-dimensional 
system with hard-wall boundary conditions in the x-direction 
and periodic boundary conditions in the y-direction. The 
solid lines for the boundary means to take hard-wall boundary 
conditions, and broken lines for the boundary means to take 
periodic boundary conditions. 



In this system the ^-component of the total momen- 
tum is not conserved, so the total number of the zero- 
Lyapunov exponents is 4. Fig. ^ is a small positive 
region of the Lyapunov spectrum normalized by the max- 
imum Lyapunov exponent Ai « 1.30, including its step- 
wise structure, while its whole positive part of the Lya- 
punov spectrum is given in the inset. The stepwise struc- 
ture of the Lyapunov spectrum is clearly different from 
the model of the previous subsection, and consists of 
two-point steps interrupted by isolated single Lyapunov 
exponents. (We call the isolated single Lyapunov ex- 
ponents interrupting the two-point steps the "one-point 
steps" in this paper from now on, partly because they are 
connected to the two-point steps in the model of Sect. 
Ill as will shown in this subsection.) As discussed in 



Sect. Ill the Lyapunov spectrum for the model with pe- 



riodic boundary conditions in both directions has two- 
point steps and four-point steps, and it is remarkable 
that to adopt the hard-wall boundary conditions in the 
x-direction and to destroy the total momentum conserva- 
tion in this direction halves the step- widths of both kinds 
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of steps. 



FIG. 17: Stepwise structure of the Lyapunov spectrum nor- 
malized by the maximum Lyapunov exponent for the quasi- 
one-dimensional system with hard-wall boundary conditions 
in the :r-direction and periodic boundary conditions in the y- 
direction. Inset: Full scale of the normalized Lyapunov spec- 
trum. The circle dots filled by black (gray) are the Lyapunov 
exponents accompanying wavelike structure (a constant be- 
havior) of the Lyapunov vector components {Sq^} partly 
shown in Fig. [18l 



FIG. 18: Time-averaged Lyapunov vector components 
{SQyj) corresponding to the n-th Lyapunov exponent A„, 
n — 149, 148, 145 and 142, as functions of the time-averaged 
particle position {q X j ) / L x normalized by the system length 
L x . The system is the quasi-one-dimensional system with 
hard-wall boundary conditions in the a;-direction and periodic 
boundary conditions in the y-direction, and the corresponding 
Lyapunov exponents are shown as the black- and gray-filled 
circles in Fig. The numerical data are fitted by a constant 
and sinusoidal functions. 



Fig. [T^ is the graphs of the time-averaged Lyapunov 
vector components {Sq^) as functions of the normal- 
ized position (q X j)/L x (the graph [5i]), corresponding 
to the one-point steps. We can clearly recognize wavclike 
structures in these graphs, and in this sense the one-point 
steps in this model should be strongly related to the two- 
point steps in the model of Sect. The Lyapunov ex- 
ponents accompanying wavelike structure of this kind of 
graphs are shown in Fig. 17 as black-filled circle dots, 
meaning that they are the one-point steps. In Fig. [l7] 
we also filled the circle dots by gray for one of the zero- 
Lyapunov exponents in which the graph of the Lyapunov 
vector components {$q^f) as a function of the position 
(q x j) is constant approximately. However it is important 
to note that the wavelength of the wave corresponding to 
the i-th one-point step in this model is 2/i, not 1/i like in 
the model of Sect. III. We fitted the graphs correspond- 
ing to the Lyapunov exponents A„, n = 149,148,145 
and 142 by the functions y — 0:149, CV148 cos(7r:r + 
Pus), «145 cos(27rx -(- (3 145 ), CV142 cos(37r:r + /?i42) respec- 
tively, with a„ and (3 n as fitting parameters. Here 
the values of the fitting parameters are chosen as 
ai49 = 0.11547, (a 14 8,/3i48) = (0.16213,-0.0089778), 
(animus) = (-0.16202,-0.01007) and (a li2 ,(3 li2 = 
(-0.15945,0.0089148). The graphs are very nicely fitted 
by a constant or the sinusoidal functions, and lead to the 
form 

fg£«W fc cos + (5) 

j = 1,2, ••• ,N, of the Lyapunov vector component 
5<£$ corresponding to the Lyapunov exponents \ n (k) 
in the fc-th one-point step with constants a' k and f3' k . 

Now we investigate the remaining steps, namely the 
two-point steps of the Lyapunov spectrum. Correspond- 
ing to these two-point steps of the Lyapunov spectrum, 
the graphs of the quantity (5q!ff /Pyj)t as functions of the 
normalized position (q X j)t/L x and the collision number 
nt (the graph [B2]) show spatial wavelike structures os- 
cillating in time. It is shown in Fig. [19 for those graphs 



corresponding to Lyapunov exponents (indicated by ar- 
rows in Fig. [I?]) in different two-point step, in the same 
collision number interval [222000,232200]. The spatial 
wavelength of the waves corresponding to the i-th two- 
point step is 2/i, which is twice as long as the wavelength 
of waves of the four-point steps in the models in Sect. 



Ill 



and the previous subsection. The period of time- 
oscillation of the wave corresponding to the z-th two- 
point step of the Lyapunov spectrum is approximately 
given by T^/i with a constant Tq. This kind of graph 
corresponding to one of the zero-Lyapunov exponents, 
namely A150, is almost constant. 



FIG. 19: Local time-averaged quantities (Sq^ /Pyj)t as 
functions of the normalized position {q X j)t/L x and the col- 
lision number nt corresponding to the Lyapunov exponents 
A n , n = 147, 144 and 141 in the same collision number inter- 
val [222000, 232200]. The system is the quasi-one-dimensional 
system with hard-wall boundary conditions in the z-direction 
and periodic boundary conditions in the y-direction, and 
the correspon ding Lyapunov exponents are indicated by ar- 
rows in Fig. O. Contour plots on the bottoms of these 
three-dimensional plots are given by red dotted lines, black 
solid lines and blue broken lines corresponding to the values 
{Sq^ /p yj )t = 0.08,0 and -0.08, respectively. 



It is important to note a relation between the time- 
oscillating period 7q of the preceding two models and 
the time-oscillating period 7q of the model in this sec- 
tion. Noting that in Figs. ||(a), |lj(a) and |l^(a) we 
plotted about one period of the time-oscillation of the 
wavelike structures of the quantities (5q^ /p y j)t in the 
collision time intervals [543000,548100], [390600,396000] 
and [222000,232200], respectively, we can get an approx- 
imate relation Tq 27o. 

The next problem is to investigate the graphs of the 
quantities (<^™ /Pyj)t as a function of the position (q x j)t 
and the collision number nt in the same two-point step 
of the Lyapunov spectrum. Fig. ^ is the contour plots 
of such graphs for the first two-point step consisting of 
the Lyapunov exponents A„, n = 147 and 146, which is 
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FIG. 20: Contour plots of the local time-averaged quan- 
tities {Sqjjjj /p y j)t as functions of the normalized position 
{<lxj)t/L x and the collision number nt corresponding to the 
first two-point step consisting of the Lyapunov exponents 
X n , n — 147 and 146 in the same collision number interval 
[212400,232200]. The system is the quasi-one-dimensional 
system with hard-wall boundary conditions in the ^-direction 
and periodic boundary conditions in the y-direction, and 
the corresponding Lyapunov exponents are indicated by the 
brace under circles in Fig. |l^. Here red dotted lines, black 
solid lines and blue broken lines correspond to the values 
{8q$/p v j)t = 0.08,0 and -0.08, respectively. 



point steps in this model correspond to the four-point 
steps in the models of the previous subsection and Sect. 
[II, except for differences in the values of their wave- 
lengths and time-oscillating periods. After all we get a 
conjecture that the Lyapunov vector components dqljj^ 

and S^yj^ ^ corresponding to the Lyapunov exponents 
constructing the fc-th two-point step are approximately 
expressed as 



+ <M,'" I /P,j>, > 
m <M,f>/P,j>, = 
— <5ry"//>„>, < 
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FIG. 21: Schematic illustration of the phase relations 
among the time-oscillating wavelike structures of the quan- 
tities (5<lyj /Pyj)t as functions of the position (q X j)t and the 
collision number n t , corresponding to the first two-point step 
in the same collision number interval. Thick gray lines means 
node lines, the region indicated by a red plus sign (+) is the 

region where the quantity (Sq^ /p y j)t is positive, and the re- 
gion indicated by a blue minus sign ( — ) is the region where 

the quantity (Sq^ /pyj)t is negative. The phase of [PI'] is 
shifted in time from the phase of [P2'\ by 7r/2. 



irk \ / 

a' k p VJ cos ( —q xj + (3' k j cos ( — n t +% 



irk 



nk 



a k p yj cos \j-q X j+f3' k j sin n t +j k j j , (6) 

j = 1, 2, • • • , N with constants a' k , f3' k and j' k , noting the 
time-oscillating period Tq w 2%. 



C. The case of hard- wall boundary conditions in 
both directions 

The last model is the case of hard- wall boundary con- 
ditions in both directions (the boundary case [-4^]). A 
schematic illustration of this system is given in Fig. |22| 
in which the solid line of the boundary means to take 
hard-wall boundary conditions. 



o 



indicated by a brace under circles in Fig. [17|, in the same 
collision number interval [212400,232200]. It should be 
noted that the positions of the nodes of two spatial 
waves belonging to the same two-point step almost co- 
incide with each other. However the phases of the time- 
oscillations of the amplitudes of the waves are shifted by 
about 7r/2 with each other. The phase relations of the 
graphs |2^(a) and (b) is visualized in the schematic illus- 
tration given in Fig. [2l| of the time-oscillating wavelike 
structures of the quantities (Sq^ IVyj)t corresponding to 
the first two-point step in the same collision number in- 
terval. Here the phase [PI '] and [P2'] correspond to Figs, 
^(a) and (b), respectively. A similar investigation of the 
Lyapunov vectors shows that the time-oscillating wave- 
like structures for the second two-point steps are like the 
phases [PI] and [P2] of Fig. [ll] except that the period Tq 
in Fig. [ll] should be replaced with the oscillating period 
Tq of this model. These results suggest that the two- 



FIG. 22: A schematic illustration of a quasi-one-dimensional 
system with hard-wall boundary conditions in both directions. 
The solid line of the boundary means to take the hard-wall 
boundary conditions. 

A small positive region of the Lyapunov spectrum nor- 
malized by the maximum Lyapunov exponent Ai ~ 1.28 
is given in Fig. |2^. The graph for the full scale of the 
positive branch of the normalized Lyapunov spectrum is 
also given in the inset of this figure. In this system the 
total momentum is not conserved any more, and the total 
number of zero- Lyapunov exponents is 2. The stepwise 
structure of the Lyapunov spectrum consists of two-point 
steps only. In this model a wavelike structure in the Lya- 
punov vector component (^Qy^) as a function of the po- 
sition (q X j) (the graph [Bl]) is not observed. 

Fig. [24| is the graphs of the quantities {Sq y ^ /p y j)t as 
functions of the normalized position (q X j)t/L x and the 
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FIG. 23: Stepwise structure of the Lyapunov spectrum nor- 
malized by the maximum Lyapunov exponent for a quasi- 
one- dimensional system with hard-wall boundary conditions 
in both directions. Inset: Full scale of the normalized Lya- 
punov spectrum. 

collision number nt (the graph [52]), corresponding to 
the Lyapunov exponents A n , n = 149,147 and 145, in 
the same collision number interval [235200, 246000]. The 
corresponding Lyapunov exponents are indicated by ar- 
rows in Fig. ^3| This figure for the two-point steps of the 
Lyapunov spectrum shows a similar wavelike structure 
to the wavelike structure of the quantities (Sq y j /Pyj}t m 
the model of the previous subsection, although one may 
think that fluctuations of these graphs in this model are 
much smaller than in the previous model. This suggest 
that the two-point steps of the Lyapunov spectrum in 
Fig. |2^ are similar to the two-point steps in the model of 
the previous subsection. The wavelength of the spatial 
waves and the the periods of the time-oscillations cor- 
responding to the i-th two-point step are 2/i and T '/i, 
respectively, and the time-oscillating period of the first 
two-point step is given approximately by the same period 
Tq (w 27o) as in the model of the previous subsection. It 
may also be noted that this kind of graph corresponding 
to the zero-Lyapunov exponent A150 is almost constant. 

FIG. 24: Local time-averaged quantities (Scffi /Pyj)t as 
functions of the normalized position (q X j)t/L x and the col- 
lision number nt corresponding to the Lyapunov exponents 
X n , n — 149, 147 and 145 in the same collision number inter- 
val [235200, 246000]. The system is the quasi-one-dimensional 
system with hard-wall boundary conditions in both directions, 
and the corresponding Lyapunov exponents are indicated by 
arrows in Fig. ^j. Contour plots on the bottoms of these 
three-dimensional plots are given by red dotted lines, black 
solid lines and blue broken lines corresponding to the values 
{Sq^f /p yj )t = 0.08,0 and -0.08, respectively. 

Fig. ^5] is the contour plots of the quantities 
(SQyj /p y j)t as functions of the normalized position 
(q X j)t/L x and the collision number n% corresponding to 
the Lyapunov exponents A„ , n — 149 and 148 in the first 
two-point step of the Lyapunov spectrum in the same col- 
lision number interval [223800, 246000]. The correspond- 
ing Lyapunov exponents are indicated by the brace under 
circles in Fig. ^3|. Similarly to the previous model, the 
nodes of two spatial waves corresponding to the same 
two-point step almost coincide with each other, and the 
phase of the time-oscillation of the wave amplitudes is 
shifted by about ir/2. This also says that the phase re- 



lations of the graphs ^5|(a) and (b) is the types of the 
phases [PI'] and [P2'] in Fig. In a similar way we 
can see that the time-oscillating wavelike structures for 
the second two-point steps of the Lyapunov spectrum for 
this model are like the phases [PI] and [P2] of Fig. |ll| 
by replacing the period Tq of Fig. [ll] with the time- 
oscillating period T Q ' w 2% of this model. These sug- 
gest an approximate expression for the Lyapunov vector 
components given by Eq. (Q), for Lyapunov vector com- 
ponents Sq^ corresponding to the Lyapunov exponents 
of the A;-th two-point steps in this model. The fact that 
the only difference between the model in this subsection 
and the model in the previous subsection is the boundary 
conditions in the y-direction suggests that the one-point 
steps of the model in the previous subsection come from 
the conservation of the y-component of the total momen- 
tum. 

FIG. 25: Contour plots of the local time-averaged quan- 
tities (5q y j /Pyj)t as functions of the normalized position 
{qxj)t/L x and the collision number n t corresponding to the 
first two-point step consisting of the Lyapunov exponents 
X„, n = 149 and 148 in the same collision number interval 
[223800,246000]. The system is the quasi-one-dimensional 
system with hard-wall boundary conditions in both directions, 
and the corresponding Lyapunov exponents are indicated by 
the brace under circles in Fig. Here, red dotted lines, 

black solid lines and blue broken lines correspond to the val- 
ues {Sqi^ /p V j)t = 0.08,0 and —0.08, respectively. 



V. CONCLUSION AND REMARKS 

In this paper we have discussed numerically the step- 
wise structure of the Lyapunov spectra and its corre- 
sponding wavelike structures for the Lyapunov vectors in 
many-hard-disk systems. We concentrated on the quasi- 
one-dimensional system whose shape is a very narrow 
rectangle that does not allow exchange of disk positions. 
In the quasi-one-dimensional system we can get a step- 
wise structure of the Lyapunov spectrum in a relatively 
small system, for example, even in a 10-particle system, 
whereas a fully two-dimensional system would require a 
much more particles. In such a system we have consid- 
ered the following two problems: [^4] How does the step- 
wise structure of the Lyapunov spectra depend on bound- 
ary conditions such as periodic boundary conditions and 
hard-wall boundary conditions? [B] How can we catego- 
rize the stepwise structure of the Lyapunov spectra us- 
ing the wavelike structure of the corresponding Lyapunov 
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TABLE I: The stepwise structures of the Lyapunov spectra and the associated wavelike structures of the Lyapunov vectors 
in the four boundary c ases consi dered in this paper. The boundary cases [Al], [A 2], [A3] and [A4] were considered in the 



sections III 



[V A , IV B| and [IV C , respectively. In this table, 5 is the number of points in the step (or the number of the zero 
Lyapunov exponents in the line specified by the label "A n = 0"), £ is the wavelength of spatial wavelike structure, and T is 
the period of time-oscillation of the wave, L x is the length of the quasi-one-dimensional rectangle, and To is constant. 
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vectors? To consider the problem [A] also means to in- 
vestigate the effects of the total momentum conservation 
law on the stepwise structure of the Lyapunov spectra. 
In this paper we considered four types of the bound- 
ary conditions; [Al] periodic boundary conditions in the 
x- and y-directions, [A 2] periodic boundary conditions 
in the ^-direction and hard-wall boundary conditions in 
the y-direction, [A3] hard-wall boundary conditions in 
the ^-direction and periodic boundary conditions in the 
y-direction, and [A 4] hard- wall boundary conditions in 
the x- and y-directions, in which we took the y-direction 
as the narrow direction of the rectangular shape of the 
system. With each boundary case we obtained differ- 
ent stepwise structures of the Lyapunov spectra. In each 
boundary condition we also considered graphs of the fol- 
lowing two quantities; [Bl] the y-component 5q y j of the 
spatial coordinate part of the Lyapunov vector of the 
j-th particle corresponding to the Lyapunov exponent 
X n as a function of the x-component q X j of the spatial 
component of the j-th particle, and [B2] the quantity 
5q y n ^ /pyj with the y-component p y j of the momentum 
coordinate of the j-th particle as a function of the po- 
sition q x j and the collision number n t . These quanti- 
ties Sq y ^ and Sq y j I Pyj give constant values in some of 
zero-Lyapunov exponents, at least approximately. We 
found that the steps of the Lyapunov spectra accompany 
a wavelike structure in the quantity Sq^ or 8q~"j /p y j, 
depending on the kind of steps of the Lyapunov spec- 
tra. A time-dependent oscillating behavior appears in 
the wavelike structure of the quantity Sq y j I Pyj, whereas 
the wavelike structure of the quantity bq^ is essentially 
stationary. Fluctuations of these quantities 5qif and 



5q y j I Pyj disturb their clear oscillatory structures, so we 
took a time-average of these quantities (a local time- 
average for the quantity Sq y " ' /p v j because of its time- 
oscillating behavior, and a longer time-average for the 
quantity Sq y n ) because it is much more stationary in time 

than the quantity 5q y j /p y j) to get their dominant wave- 
like structures. In Table. | we summarize our results 
about the characteristics of the stepwise structures of 
the Lyapunov spectra and the wavelike structures of the 
Lyapunov vectors in the four boundary cases [Al], [A 2], 
[A3] and [A4]- In Lyapunov exponents in each step of the 
Lyapunov spectra the wavelike structures of the quantity 



^Qyj* 01 ^Qyj' /Pyj are approximately orthogonal to each 
others in space (in the sense of Eq. @), in space and 
time (in the sense of Eq. (Q)) or in time (in the sense of 
Eq. (^|)), and this fact suggests that the wavelike struc- 
ture of these quantities are sufficient to categorize the 
stepwise structure of the Lyapunov spectra in the quasi- 
one-dimensional systems considered in this paper. 

Different from a purely two-dimensional model such as 
a square system in which each particle can collide with 
any other particle, in the quasi-one-dimensional model a 
separation between the stepwise region and the smoothly 
changing region of the Lyapunov spectrum is not clear. 



» 



In Rcf. 1 17 1 this point was explained as being caused by 
the fact that particles interact only with the two nearest- 
neighbor particles, whereas in the purely two-dimensional 
low-density system particles can interact with more than 
two particles. 

In this paper we considered the quasi-one-dimensional 
systems only. However there should be many interesting 
other situations in which we can investigate structures of 
the Lyapunov spectrum and the Lyapunov vectors. For 
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example, we may investigate the effect of the rotational 
invariance of the system on such structures by consider- 
ing a two-dimensional system with a circle shape. One 
might also investigate the system in which the orbit is 
not deterministic any more, in order to know whether the 
deterministic orbit plays an important role in the step- 
wise structure of the Lyapunov spectrum or not. It may 
also be important to investigate the dependence of the 
stepwise structures of the Lyapunov spectra on the spa- 
tial dimension of the system, for example, to investigate 
any structure of the Lyapunov spectra for purely onc- 
or three-dimensional many particle systems. (Note that 
the quasi-onc-dimensional systems considered in this pa- 
per are still two-dimensional systems in the sense of the 



phase space dimension.) These problems remain to be 
investigated in the future. 
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